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The random matrix theory is used to bridge the network structures and the dy- 
namical processes defined on them. We propose a possible dynamical mechanism for 
the enhancement effect of network structures on synchronization processes, based 
upon which a dynamic-based index of the synchronizability is introduced in the 
present paper. 
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The impact of network structures on the synchronizability of the identical 
oscillators defined on them is an important topic both for theory and potential 
applications. From the view point of collective motions, the synchronization 
state is a special elastic wave occurring on the network, while the initial state 
is a abruptly assigned elastic state. The synchronizability should be the tran- 
sition probability between the two states. By means of the analogy between 
the collective state and the motion of an electron walking on the network, we 
can use the quantum motion of the electron to find the motion characteristics 
of the collective states. The random matrix theory (RMT) tells us that the 
nearest neighbor level spacing distribution of the quantum system can capture 
the dynamical behaviors of the quantum system and the corresponding classical 
system. A Poison distribution shows that the transition can occur only between 
successive eigenstates, while a Wigner distribution shows that the transition 
can occur between any two eigenstates. A Brody distribution, an intermediate 
between the two extreme conditions, can give us a quantitative description of 
the transition probability. Hence, it can be used as an index to represent the 
synchronizability. As examples, the Watts-Strogatz (WS) small-world networks 
and the Barabasi- Albert (B A) scale-free networks are considered in this paper. 
Comparison with the widely used eigenratio index shows that this index can 
describe the synchronizability very well. It is a dynamic-based index and can 
be employed as a measure of the structures of complex networks. 



I. INTRODUCTION 



Recent years witness an avalanche investigation of complex networks |l|,|2j,y]. Complex 
systems in diverse fields can be described with networks, the elements as nodes and the 
relations between these elements as edges. The structure-induced features of dynamical 
systems on networks jattract special attentions, to cite examples, the synchronization of 
coupled oscillators 
to external stimuli lUi 




, the epidemic spreading [9 



10| and the response of networks 
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Synchronization is a wide-ranging phenomenon which can be found in social, physical 
and biological systems. Recent works show that some structure features of complex net- 
works, such as the small-world effect and the scale-free property, can enhance effectively the 
synchronizabilities of identical oscillators on the networks, i.e., synchronization can occur in 
a much more wide range of the coupling strength. 

We consider a network of N coupled identical oscillators The network structure can 
be represented with the adjacent matrix A, whose element Ajj is and 1 if the nodes i and 
j are disconnected and connected, respectively. Denoting the state of the oscillator on the 

node i as x l = I ' I ,the dynamical process of the system is governed by the following 



Pi 



equations, 



N 

± l = F{x i )-aY,L l3 Q{x') ) (1) 

3=1 

where x l = F(x l ) governs the individual motion of the ith oscillator, a the coupling strength 
and Q{xi) the output function. The matrix L is a Laplacian matrix, which reads, 



Lij 



N 

A%s = ki (i = j / 



(2) 



-Ai (i j) 

ki&ij A%j 5 

where ki is the degree of the node i, i.e., the number of the nodes connecting directly with 
the node i. The eigenvalues of L are real and nonnegative and the smallest one is zero. 
That is, we can rank all the possible eigenvalues of this matrix as = 71 < 72 < ■ • • < Jn- 
Herein, we consider the fully synchronized state, i.e., \x l (t) — x^(t)\ — > as t — ► 00 for any 
pair of nodes i and j. 

Synchronizability of the considered network of oscillators can be quantified through the 
eigenvalue spectrum of the Laplacian matrix L . Here we review briefly the general frame- 



work established in 



kJ. The linear stability of the synchronized state is determined 



by the corresponding variational equations, the diagonalized iV block form of which reads, 
z = [DF(s) + 'yDQ(s)] z. z is the different modes of perturbation from the synchronized 
state. For the ith block, we have 7 = a~fi,i = 1, 2, • • • ,N. The synchronized state is stable if 
the Lyapunov exponents for these equations satisfy T(a'ji) < for i = 2, 3, • ■ ■ ,N. Detailed 



investigations 



12. 13] show that for many dynamical systems, there is a single interval of 



the coupling strength (01,02), in which all the Lyapunov exponents are negative. In this 
case, the synchronized state is linearly stable if and only if ^ < ^ = x- While x depends 
on the the dynamics, the eigenratio R = depends only on the topological structure of 
the network. Hence, this eigenratio represents the impacts of the network structure on the 
networks's synchronizability. This framework has stimulated an avalanche investigation on 
the synchronization processes on complex networks. It has been widely accepted as the 
quantity index of the synchronizability of networks. 

However, the eigenratio is a Lyapunov exponent-based index. It can guarantee the linear 
stability of the synchronized state. It can not provide enough information on how the 
network structure impacts the dynamical process from an arbitrary initial state to the final 
synchronized state. How the structures of complex networks impact the synchronization is 
still a basic problem to be understood in detail. In this paper, by means of the random matrix 
theory (RMT) , we try to present a possible dynamical mechanism of the enhancement 
effect, based upon which we suggested a new dynamic-based index of the synchronizabilities 
of networks. 



II. DYNAMIC-BASED INDEX OF SYNACHRONIZ ABILITY 



The RMT was developed by Wigner, Dyson, Mehta, and others to understand the energy 



levels of complex quantum systems, especially heavy nuclei [14j . Because of the complexity of 
the interactions, we can postulate that the elements of the Hamiltonian describing a heavy 
nucleus are random variables drawn from a probability distribution and these elements 
are independent with each other. A series of remarkable predictions are found to be in 
agreement with the experimental data. The great successes of RMT in analyzing complex 
nuclear spectra has stimulated a widely extension of this theory to several other fields, such 



as the quantum chaos, the time series analysis 15^ 



mesoscopic systems, the complex networks 



20, 



16 



21. 



17L 

10 



18, 



1£ [jthe transport in disordered 



25 



27 



2a ]. and even the 



QCD in field theory. For the complex quantum systems, the predictions represent an average 
over all possible interactions. The deviations from the universal predictions are the clues 
that can be used to identify system specific, non-random properties of the system under 
consideration. 
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One of the most important concepts in RMT is the nearest neighbor level spacing (NNLS) 
distribution 14] . Enormous experimental and numerical evidence tells us that if the classical 
motion of a dynamical system is regular, the NNLS distribution of the corresponding quan- 
tum system behaves according to a Poisson distribution. If the classical motion is chaotic, 
the NNLS distribution will behave in accordance with the Wigner-Dyson ensembles, i.e, 
~ s -exp(— us 2 ), s is the NNLS. The NNLS distribution of a quantum system can tell us the 
dynamical properties of the corresponding classical system. This fact is used in this paper 
to bridge the structure of a network with the dynamical characteristics of the dynamical 
system defined on it. 

From the state of the considered system, X = j J \z = 1, 2, • • • , N , we can construct 

w 

the collective motion of the system as, 



$(M) = x i (9 i ,I i ,t,Y),i= 1,2,- ■■ ,N, (3) 

where Qi and Ii are the phase and the amplitude of the oscillator i. Y is the other oscillation- 
related parameters. $ describes the elastic wave on the considered network and Eq(3) 
presents the displacements at the positions i = 1, 2, • • • , N at time t. Because of the iden- 
tification of the oscillators, the individual motions should behave same except the phases 
and the amplitudes. The synchronization process can be described as the transition from 

( Qi o°2 ■ ■ ■ o% , 

an arbitrary initial collective state, , to the final fully synchronized state, 

OT n \ _ ( 9 s 2 yn \ ( 6 S j» n 

jsyn J ~ yjsyn J " " ^ jsyn 

The probability of the transition should be the synchronizability of the considered network. 
The larger the transition probability, the easier for the system to achieve the fully synchro- 
nized state. 

The collective states are the elastic waves on the considered network. This kind of 



electron wal 


ring 


29 


U 


31 







this analogy is used to extend the concept of Anderson localization state to the classical 
phenomena as elastic and optical waves. In this paper we will use it to find a quantitative 
description of the transition probability between the collective states. 
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The tight-binding Hamiltonian of an electron walking on the network reads, 

N N 

H = Y,e i -\i){i\+Y^A ij t ir \i){j\, (4) 

i=l i^j 

where is the site energy of the zth oscillator, the hopping integral between the nodes i 
and j. Because of the identification of the oscillators, all the site energies are same, denoted 
with Si = e. Generally, we can set e = and Uj = 1, which leads to the relation H — A. 
Ranking the spectrum of A as Ai > A2 > ■ • ■ > \n, we denote the corresponding quantum 
states with \Pj \i = 1, 2, 3, • • • , N . Hence, the NNLS distribution of the adjacent matrix A 
can show us the dynamical characteristics of the collective motions. 

If the NNLS obeys the Poisson form, the transition probability between two eigenstates 
and \P n will decrease rapidly with the increase of |A m — A n |, and the transition occurs 
mainly between the nearest neighboring eigenstates. This state is called quantum regular 
state. If the NNLS obeys Wigner form, the transitions between all the states in the same 
chaotic regime the initial state belongs to can occur with almost same probabilities. The 
electron is in a quantum chaotic state. 

The corresponding collective states of the classical dynamical system to the quantum 
chaotic and regular sates are called collective chaotic and collective regular states, respec- 
tively. If the dynamical system is in a collective chaotic state, the collective motion modes in 
same chaotic regimes can transition between each other abruptly, while if the system is in a 
collective regular state only the neighboring collective motion modes can transition between 
each other. Generally, a dynamical system may be in an intermediate state between the 
regular and the chaotic states, which is called soft chaotic state. 

The NNLS distribution can be obtained by means of a standard procedure. The first 
step is the so-called unfolding. In the theoretical predictions for the NNLS, the spacings 
are expressed in units of average eigenvalue spacing. Generally, the average eigenvalue 
spacing changes from one part of the eigenvalue spectrum to the next. We must convert the 
original eigenvalues to new variables, called unfolded eigenvalues, to ensure that the spacings 
between adjacent eigenvalues are expressed in units of local mean eigenvalue spacing, and 
thus facilitates comparison with analytical results. Define the cumulative density function as, 
G(X m ) = N f* m g(X)dX, where g(X) is the density of the original spectrum. Dividing G(X) 
into the smooth term G av (X) and the fluctuation term G/(A), i.e., G(X) = G av (X) + Gf(X), 
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the unfolded energy levels can be obtained as, 

If the system is in a soft chaotic state, the NNLS distribution can be described with the 
Brody form |22(], which reads, 



PU) = t ■ gP- 1 ■ exp 
V 



(6) 

,v; 

We can define the accumulative probability distribution as, Q(s) = f* P(u)du. The pa- 
rameter (3 can be obtained from the linear relation as follows, 



InT(s) = In 



In 



1 



j3lns — (3lnn. (7) 



l-Q(s) 

For the special condition (3 — 1, the probability distribution function (PDF) P(s) degener- 
ates to the Poisson form and the system is in a regular state. For another condition (3 = 2, 
the PDF obeys the Wigner-Dyson distribution P(s) oc s ■ exp(— 7 • s 2 ) and the system is in a 
hard chaotic state. If the system is in an intermediate soft chaotic state, we have, 1 < (3 < 2. 

Hence, from the perspective of random matrix theory, the synchronizability can be de- 
scribed with the parameter (3. The larger the value of (3, the easier for the system to 
become fully synchronized. By this way we find a possible dynamical mechanism for the 
enhancement effects of the network structures on the synchronization processes. 

III. RESULTS 

In reference j^, the authors prove that the spectra of the Erdos-Renyi, the Watts- 
Strogatz(WS) small-world, and the growing random networks (GRN) can be described in a 
unified way with the Brody distribution. Herein, we are interested in the relation between 
the parameter (3 and the eigenratio R. Detailed works show that R is a good measure of the 
ity of complex networks, especially the small world and scale-free networks 



syncnromza 
0,0,QI 



synchronizabi 

Figure 1 shows the relation between (3 and R for WS small- world networks j38|. We 
use the one-dimensional regular lattice-based model. In the regular lattice each node is 
connected with its d right-handed neighbors. Connecting the starting and the end of the 



8 



lattice, with the rewiring probability p r rewire the end of each edge to a randomly selected 
node. In this rewiring procedure self-edge and double edges are forbidden. Numerical results 
for WS small-world networks with N = 3000 and d = 2 are presented. We can find that the 
Brody distribution can capture the characteristics of the PDFs of the NNLS very well, as 
shown in the panel (a) in Fig.l. With the increase of p r , the parameter 3 increases rapidly 
from 1.32 to ~ 1.85, while the parameter R decreases rapidly from 3000 to ~ 10 1 . Hence, 
there exists a monotonous relation between the two parameters 3 and R. 

n 

Figure 2 gives the results for Barabasi- Albert (BA) scale-free networks [39j. Starting 
from a seed of several connected nodes, at each time step connect a new node to the existing 
graph with w edges. The preferential probability to create an edge between the new node 
and an existing node / is proportional to its degree, i.e., pba oc k(f). Numerical results 
for BA scale-free networks with N = 3000 and w = 1,2,3, •• ■ ,10 are presented. All the 
PDFs of the NNLS obey the Brody distribution almost exactly. With the increase of w, the 
parameter 3 increases from 0.35 to ~ 1.70, while the parameter R decreases from 51509 to 
~ 10 1 . We can find also a monotonous relation between the two parameters 3 and R. 

For w = 1, we have 3 = 0.35 < 1. That is, rather than the "repulsions" or un-correlations 
between the levels, there are a certain "attractiveness" between the levels. In the construc- 
tion of the BA networks with w — 1, each time only one node is added to the existing 
network. The resulting network is a tree-like structure without loops at all. Dividing the 
network into subnetworks, we can find that many of them have similar structures, which 
leads their corresponding level-structures being almost same. Because of the weak coupling 
between the subnetworks, the total level structure can be produced just by put all the corre- 
sponding levels together. This kind of level-structure will lead many NNLS tending to zero. 
Hence, 3 < 1 is an extreme case induced by tree-like structure. This special kind of tree-like 
BA networks can not enhance the synchronization at all. 



IV. DISCUSSIONS 



In summary, by means of the NNLS distribution we consider the collective dynamics in 
the networks of coupling identical oscillators. For the two kinds of networks, we can find the 
monotonous relation between the two parameters 3 and R. This monotonous relation tells 
us that the high synchronizability is accompanied with a high extent of collective chaos. The 
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collective chaos may increase significantly the transition probability of the initial random 
state to the final synchronized state. The collective chaotic processes may be the dynamical 
mechanism for the enhancement impacts of network structures on the synchronizabilities. 

The parameter (3 in the NNLS distribution can be a much more informative measure 
of the synchronizability of complex networks. It reveals the information of the dynamical 
processes from an arbitrary initial state to the final synchronized state. It can be regarded in 
a certain degree as the bridge between the structures and the dynamics of complex networks. 

One paradox may be raised about the argument in the present paper. The Wigner 
distribution implies a larger correlation between the eigenstates of the network than does 
the Poisson distribution. At the same time, one can reverse the argument that Wigner 
distribution implies level repulsion and, therefore, different frequencies of oscillation of the 
normal modes, and therefore no synchronization when these modes are coupled. It should 
be emphasized that the eigenratio R and the index (3 should be used together to capture the 
impacts of the network structures on the synchronization processes. R represents the linear 
stability of the synchronized state, but it can not tell us how the final synchronized state is 
reached from the initial state. On the other hand, (3 provides us a possible mechanism for 
this dynamical processes, but it can not tell us the transition orientation. R and (3 reflect 
some features of the impacts of the network structures on the synchronization processes, but 
there may be some new important features to be found. 
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FIG. 1: The relation of (3 versus R for the constructed WS small-world networks, (a) Several 
typical results for PDF of the NNLS. In the interested regions a brody distribution can capture 
the characteristics very well, (b) With the increase of the rewiring probability p r the eigenratio R 
decreases rapidly, (c) With the increase of the rewiring probability p r the parameter (3 increases 
rapidly, (d) The monotonous relation between the two parameters (3 and R. 
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FIG. 2: The relation of (3 versus R for the constructed BA scale-free networks, (a) Results for PDF 
of the NNLS. In the interested regions a brody distribution can capture the characteristics very 
well, (b) With the increase of w the eigenratio R decreases significantly, (c) With the increase of w 
the parameter (3 increases significantly, (d) The monotonous relation between the two parameters 
(3 and R. 



